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Phenomenological Ginzburg-Landau theory is used to calculate the possible spontaneous vortex 
states that may exist at corner junctions of d^2_y2 + ix-wme, (where x = s or x — d^y) and s- 
wave superconductors. We study the magnetic flux and the critical current modulation with the 
junction orientation angle 9, the magnitude of the order parameter, and the magnetic field H . It is 
seen that the critical current Jc versus the magnetic flux $ relation is symmetric / asymmetric for 
X = dxy/s when the orientation is exactly such that the lobes of the dominant d^2_y2-wave order 
parameter points towards the two junctions, which are at right angles for the corner junction. The 
conclusion is that a measurement of the Ic{^) relation may distinguish which symmetry {d^-2_y2 +is 
or d^2_y2 + idxy) the order parameter has. 



I. INTRODUCTION 

One of the main questions in the research activity on 
high- Tc superconductors nowadays is the identification of 
the order parameter symmetry and its underlying mech- 
anism The most possible scenario is that the bulk 
pairing state has a d^2_y2-wave character. Theoretical 
calculations, suggest that an imaginary s-wave compo- 
nent which breaks the time reversal symmetry is induced 
in some cases, wherever the d2,2_j^2-wave order parameter 
varies spatially such as near a vortex, or near the surface 
Also the observation of fractional vortices on a tri- 
angular grain boundary in YBa2Cu307 by Kirtley et al. 
[Ql , may indicate a possible violation of the time- reversal 
symmetry near grain boundary. Theoretical explanation 
of this experiment is given by Bailey et. al. in Ref. ||] 
where they study a triangular grain boundary in d-wave 
superconductors. They conclude that under the assump- 
tion of d-wave symmetry, the flux at the edges of this 
triangle can take the values ±$o/2, which does not agree 
with the experiment. However under the assumption of 
dr^2_y2 -|- is- wave symmetry an intrinsic phase shift (j)c{x) 
exists in each triangle edge. In turn the phase (l){x) must 
change in order to connect the different values of 0c in 
each segment. This arrangement leads to fractional vor- 
tices or antivortices at each three corners, in agreement 
with the experiment. 

Another pairing state which breaks the time reversal 
symmetry is the dx^_y2 -\- id^y-'wave. Patches of complex 
dxy components are induced around magnetic impuri- 
ties at low temperatures in a dj.2_j^2-wave superconductor 
forming a phase coherent state as a result of tunneling 
between different patches Violation of parity and 
time reversal symmetry occurs in this state. Also on the 
high field region, H < Hc2 the d^2_y2-wave state can be 
perturbed by the external filed, producing a dx2-y2 + idxy 
state in the bulk @] . 



The observation of the splitting of the zero energy peak 
in the conductance spectra at low temperatures indicates 
that a secondary component is induced which violates lo- 
cally the time reversal symmetry Theoretical expla- 
nation based on surface- induced Andreev states, has been 
proposed ^ . Recently the field dependence of this split- 
ting has been observed in the tunneling spectra of YBCO 
|lO| , p| . This observation is consistent with a dj;2_y2 +is 
surface order parameter as well as a dx2-y2 + idxy bulk 
order parameter. Another question which can be asked 
is to what extend, the observation of a symmetric mag- 
netic interference pattern in the corner junction experi- 
ments is an identification of (i2,2„j,2-wave symmetry, 
or could also imply a dx^-y^ + idxy pairing state also? 
In this work we propose a phase sensitive experiment 
based on the Josephson effect, which may be used to dis- 
tinguish which symmetry [dx^-y'i + is or dx2-y2 + idxy) 
the order parameter has near the surface. We study the 
static properties of a frustrated junction which is made 
of two one-dimensional junctions, of dx2-y2 + ix-wave, 
(where x = dxy or x ~ s) and s-wave superconductors. 
By introducing an extra relative phase in one part of 
this junction, the above junction can be mapped into 
the corner junctions experiments |P,p^. We examine the 
spontaneous flux and the critical current modulation of 
the vortex states with the junction orientation angle 9, 
the magnitude of the secondary component Ug, and the 
magnetic field H. In each case we derive simple argu- 
ments which are useful to discriminate between the time 
reversal symmetry broken states. For example, when the 
orientation is exactly such that the lobes of the dominant 
d2.2_j,2-wave order parameter points towards the junction 
interface the magnetic interference pattern is symmet- 
ric (asymmetric) when the secondary order parameter is 
X = dxy{s). This is verified for small junctions as well as 
in the long junction limit, and can be used to distinguish 
between broken time reversal symmetry states. 
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FIG. 1. (a) A single Josephson junction between super- 
conductors A and B with a two component order parameter. 
The angle between the crystalline a axis of A and the junc- 
tion interface is 6. (b) The geometry of the corner junction 
between a mixed symmetry superconductor, and an s-wave 
superconductor. 
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Here 0^^^'' is the phase of the order parameter in su- 
perconductor A{B). Then phcnomcnological Ginzburg- 
Landau theory is used to calculate the supercurrent den- 
sity given by 1 13 
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m* ,1711,1711 are the effective masses that enter into the 
Ginzburg-Landau equations. In the following these 
masses are taken equal to an effective mass to*. 

We restrict to the case where B is s-wave. In this case 
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0, and =constant. We define 
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the relative phase difference between the two supercon- 
ductors. We consider the case where the intrinsic phase 
difference within superconductor A is 4>2 — 4>i = 12. 
Then the order parameter in A is complex and breaks 
the time reversal symmetry. The supercurrent density 
can be written as: 



The rest of the paper is organized as follows. In Sec. II 
we discuss the Josephson effect between a superconduc- 
tor with broken time reversal symmetry and an s-wave 
superconductor. In Sec. Ill the geometry of the corner 
junction is discussed. In Sec. IV we present the results 
for the magnetic flux of the spontaneous vortex states in 
corner junctions with some intrinsic magnetic flux. In 
Sec. V the parameters which can modulate the sponta- 
neous flux and the critical currents are considered. In 
Sec. VI a connection with the experiment is made. Fi- 
nally, a summary and discussion are presented in the last 
section. 



II. JOSEPHSON EFFECT BETWEEN TWO 
SUPERCONDUCTORS WITH MIXED WAVE 
SYMMETRY 

We consider the junction shown in Fig. 0(a), where 
two superconductors {A in the region z > t and B in the 
region z < Q), are separated by an intermediate layer. We 
assume that each superconductor has a two component 
order parameter. The order parameter for each compo- 
nent k[k — 1,2) in the superconductors, can be written 
as 



Jict)) = JcSm{(j> + (l)c), 
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where Ji = Jc2i-, J2 — ^Jc22- The Josephson critical 
current density Jc is scaled in units of J^o — i^^- Two 
special cases are the following: 

i) For d,j.2_y-2 + is wave case the magnitude of the 
dr^2_y2-wave component in (Q) is = nio cos(20) , where 
6 is the angle of the crystalline a-axis of superconduc- 
tor A with the junction interface. The magnitude of 
the secondary order parameter in superconductor A is 

= «20 = O.lniQ. 

ii) For d^2_y2 + idxy wave case, the magnitude of 
the dj.2_j,2-wave component in (Q) is given by — 
niocos(20), while the dxy wave component is = 
n2osin(20), where 77.20 = O.lnio- This order parameter 
can occur in the following way: The order parameter 
magnitude for the d-wave state Ao(6') = Aq cos(20) is 
an equal admixture of pairs with orbital moment = 



2 



±2, and can be written as Ao{6) = (Ao/2)[exp(2i6') + 
exp(— 2i0)]. In the presence of perturbation such as 
(ferromagnetically) ordered impurity spins Sz the coef- 
ficients of Lz = ±2 components will shift linearly in Sz 
with opposite signs. The final state will be Ao{9) — > 
Ao{0) + iSzAi{9), where Ai{0) = sm(20). The strength 
of the secondary component is proportional to the per- 
turbation Sz- 



III. THE CORNER JUNCTION GEOMETRY 



We consider the corner junction shown in Fig. |1|( 
between a superconductor with broken time reversal sym- 
metry at the surface and an s-wave superconductor. If 
the angle of a-axis with the interface in the a;-dircction is 
9, then the corresponding angle in the ^-direction will be 
tt/2 — 6. We map the two segments each of length L/2 
where L = lOXj of this junction into a one-dimensional 
axis. In this case the two dimensional junction can be 
considered as being made of two one dimensional junc- 
tions described in Sec. II connected in parallel. Their 
characteristic phases (j)ci and (I>c2 depend upon the angle 
9. We call this junction frustrated since the two segments 
have different characteristic phases 0ci, </'c2- The fabrica- 
tion details of corner junctions or superconducting quan- 
tum interference device (SQUID), between sample faces 
at different angles can be found in Ref. PJl^]. 

The superconducting phase difference across the 
junction is then the solution of the sine-Gordon (s-G) 
equation 



JcSin[0(x) +</)e(x)] 



with the boundary conditions 



dx 
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The length x is scaled in units of the the Josephson pen- 
etration depth given by 



hc^ 



inedJco 



where d is the sum of the s-wave, and mixed wave Xab 
penetration depths plus the thickness of the insulator 
layer. The relative phase (j)c{x) is (f>ci{4'c2) in the left 
(right) part of the junction. The external magnetic field 
H, scaled in units of He — 2edx, applied in the y di- 
rection, which is considered small compared to Aj. The 
bias current per unit length 1°^ in the overlap geometry 
is scaled in units of -^Hc, and is uniformly distributed 
along the entire x axis of the junction. 

We can classify the different solutions obtained from 
Eq. (0) with their magnetic flux content 
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stable solutions (j> = — 0ci -I- 2ni7r 
(0 = — 0C2 + 2n27r), for Ui = 0,-1, i = 1,2, that exist in 
the left(right) junction, of d^2_y2 + is- wave and s-wave su- 
perconductors, when considered uncoupled, at zero current, 
versus the orientation angle 6. Each junction has length L/2, 
where L = lOAj, and 0ci, {4>c2) is the extra phase difference 
in the left (right) junction due to the different orientations. 
The arrows denote the variation of the phase (j) in order to 
connect these stable solutions in the frustrated junction ge- 
ometry. We present three possible solutions i.e. n — 0, —1, 1, 
and down(up) arrow denotes negative(positive) magnetic flux. 



1 

2^< 



(9) 



where (pi^^L) is the value of <j) at the right (left) edge of 
the junction, in units of the flux quantum = • 



IV. SPONTANEOUS VORTEX STATES 

Firstly let us consider the case where the two one- 
dimensional junctions of dx2_y2 + ix-wave where x — s 
OT X = dxy, and s-wave supeconductors, each of length 
L/2, described in Sec. II are uncoupled. Then for 
< X < L/2 the stable solutions for the s-G equation 
are 4>{x) — — 0ci + 2ni7r, where ni = 0, ±1, ±2, while 
for i/2 < a; < L the stable solutions for the s-G equation 
are (f>{x) = —(j)c2 + 2n27r, where n2 — 0, ±1, ±2, where 
(f'ci, 4'c2, are the relative phases in each part of the junc- 
tion due to different orientations. These solutions are 
plotted in Fig. ^, for n,; = 0, — 1, i = 1, 2 as a function of 
the orientation angle 9. When the frustrated junction is 
formed, and we consider the above junctions in parallel, 
the phase (j) is forced to change around x = L/2, to con- 
nect these stable solutions. This variation of the phase 
(j), along the junction describes the Josephson vortices. 
The flux content of these states (in units of <I>o) is ||lj] 

$ = [(t>{L) ~ (t>{Q)]/2n = (-0,2 + (t>ci + 2mi)/2Ti, 

(10) 

where the n-value (n = ni — n2 = 0, ±1,±2,...) dis- 
tinguishes between solutions with different flux content. 
We will concentrate to solutions called modes with the 
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TABLE I. The magnetic flux ($) in terms of (j)ci, (j>c2 for 
the spontaneous solutions that exist in the corner junction 
geometry between a superconductor with time reversal broken 
symmetry and an s-wave superconductor {<j)c\ , 4>c2 is the extra 
phase difference in the two edges of the corner junction due 
to the different orientations, of the a-axis of the dominant 
dj.2_y2-wave superconductor). We present only the minimum 
flux states n = 0, —1, 1. 



Vortex state n 


Magnetic flux ($) 





{~4'c2 + <^cl)/27r 
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(-0C2 + (t>ci + 27r)/27r 


-1 


{-(I>c2 + (t>ci - 2n)/2-K 



minimum flux content i.e., n — 0, 1, —1. Their magnetic 
flux in terms of (pel , 4'c2 is shown in table |. Generally 
the flux content is fractional i.e. is neither integer nor 
half-integer, as a consequence of the broken time reversal 
symmetry of the problem. 

In the actual numerical simulations, the stable solu- 
tions of the sine-Gordon equation in the left (right) part 
of the junction are taken as the initial conditions for the 
iteration procedure. For example for the n — Q solu- 
tion the phase 0(x) is taken 0(x) = — 0ci {—4>c2) in the 
left (right) part of the junction, as an initial condition 
and then is iterated until convergence. Besides if we take 
as initial condition , 4){x) — ~(f>ci, in the left side, and 
(j){x) = —2tt — (j)c2 in the right side, the final state of 
the system, after the iteration procedure, is the solution 
which we call n = —1, with negative magnetic flux, and 
not exactly opposite to n = 0. We comment here that 
the solutions after the iteration procedure have smooth 
variation as a function of the position, as opposed to the 
step function variation of the initial conditions. 

For the — junction (pd = 4>c2 = 0, and the flux 
becomes $ = n, so we say that the flux is quantized in 
integer units of $0- In this case, there exist solutions with 
flux $ = -1,0, 1, ... These solutions, when n^O 
are stabilized by the application of an external magnetic 
field. In the case of a junction with some spontaneous 
fiux, at least for the modes with lower flux content, the 
external field is not necessary since the spontaneous mag- 
netization state is stable. 

In the case of — tt junction, where the intrinsic phase 
in the right (left) part of the junction is (/ic2 = — (0c1 ~ 
0), the stable solutions of the s-G equation are (j){x) — 
2mi for the left part, while 4>{x) = 7r(2n -I- 1) for the 
right part of the junction. In this case a — tt junction is 
formed. The corresponding flux becomes $ = (n+l/2)7r, 
and the particular values ofn = 0,n==— 1 give the 
half vortex and antivortex solutions, with opposite fluxon 
content, <& = 0.5 and $ = —0.5 respectively. 



V. MAGNETIC FLUX AND CRITICAL 
CURRENT MODULATION 

In the following we will describe three parameters 
which can alter the spontaneous flux and the critical 
currents of the vortex states described in the previous 
section, in a corner junction between a superconductor 
with time reversal broken symmetry and an s-wave su- 
perconductor. These include the orientation angle 9, the 
magnitude of the secondary order parameter Us , and the 
magnetic field H. In each parameter separately we will 
point out the differences between the dx2_y2 + is-wave, 
and dx2_y2 -f idxy-'^Sive. 

A. Junction orientation 

For the dx2_y2 -f is-wave case, we consider first the sit- 
uation where 9 is varied from to 7r/2. In Fig. || we plot 
the spontaneous magnetic fiux versus the angle {9) for 
the different modes n ~ 0, —1, 1 in the corner junction 
geometry. As we can see the magnetic flux changes with 
orientation. For angle 9 close to or 7r/2 the spontaneous 
modes existing at _ff = are separated by an integer 
value of the magnetic flux. This is also the case in the 
pure s-wave superconductor junction problem. The dif- 
ference is that the modes are found displaced to fractional 
values of magnetic flux, contrary to the s-wave case where 
the magnetic flux takes on integer values aX H — 0. In 
particular the vortex solution in the n = mode (solid 
line) contains less that half a fluxon for 6* = 0, and as 
we increase the angle 9 towards 7r/4 it continuously re- 
duces its flux, i.e. it becomes flat exactly at 9 — 7r/4 
and then it reverses its sign and becomes an antivortex 
with exactly opposite flux content at 9 = jl from that 
at = 0. In addition we have plotted in Fig. ^ the 
phase distributions for the mode n = in different orien- 
tations = 0, 7r/4, 7r/2. The transition form the vortex to 
the antivortex mode as the orientation changes is clearly 
seen in this figure. Note that the solutions in this mode 
remain stable for all junction orientations. This is seen 
in Fig. ^ where we plot the lowest eigenvalue (Ai) of the 
linearized eigenvalue problem as a function of the angle 9 
|l5| . We see that Ai > 0, denoting stability for all values 
of the angle 9 in this mode. 

Let as now examine the solution in the n = —\ mode, 
(dotted line in Fig. We see that at = it has 

negative fiux, which in absolute value is more than $0/2 
and as we increase the angle 9 it decreases its fiux to 
a full antifiuxon when the orientation is slightly greater 
than 7r/4 and than to flux greater than $0 when 9 reaches 
7r/2. As seen in Fig. ^ this solution becomes unstable 
at a point to the left of 6* = 7r/4 (point t) due to the 
abrupt change of flux at this angle. More strictly the 
instability sets in due to the competition between the 
slope of the phase at the edges of the junction and at the 
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FIG. 3. The magnetic flux $ as a function of the angle 9, 
for the various vortex states, n = 0, —1, 1, that exist sponta- 
neously in a corner junction between a d^2_y2 + is-wave and 
an s-wave superconductor, with length L — lOAj. The flux 
for ^ = is fractional. 

junction center as the angle 9 approaches the value 7r/4. 
At this point the slope competition makes the antivortex 
unstable. This is seen in Fig. ^) (dotted line) where 
the phase distribution for the n — ~1 mode solution 
is plotted at the point where the instability starts i.e. 
e = 0.2427r. 

Finally the solution in the n — 1 mode contains more 
than one fluxon at 6* = and is clearly unstable. It be- 
comes stable at an angle slightly on the right of 9 = tt/4, 
(point v in Fig. ^) where the flux varies more smoothly, 
see 9 = 0.2587r in Fig. |4|3. At = 7r/2 it contains more 
than <I>o /2 in flux. We expect a time reversal broken sym- 
metry state like d^^-y^ +is to be characterized by either 
the solution in the fractional vortex or antivortex mode, 
because due to the different character of these solutions 
a change from one variant to the other would demand 
the application of an external current or magnetic field 
and in this sense it would cost additional energy. So since 
these states are stable in external perturbations, once the 
system is prepared in one of these it will remain to that 
state. 

In general we see that for each value of 9 there exist in 
the junction a pair of stable solutions which when apply- 
ing an external bias current will lead to observable critical 
currents. In Fig. ^ we plot the overlap critical current 
per unit length /™ as a function of 9, a,t H — 0, for the 
n = 0,-1, 1-mode solutions, in the (1^2 _y2 + is-waye case. 
In the overlap geometry the current is distributed in the 
entire x-axis. In the calculations we have taken into ac- 
count that the Josephson critical current density Jc has 
a characteristic variation with the orientation. We find 
that for a given orientation it is possible for the junction 
current density to vary in the way that several modes 
with different critical currents can exist. In Fig. ^ we 
plot the current density when the total current is max- 
imum, for different modes, and orientations, which will 
give us information about the actual shapes of the vor- 
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FIG. 4. The phase distribution of the vortex solutions a) 
n = 0, at e = 0, 7r/4, 7r/2; b) n = -1, at 6* = 0, 0.2A2n, where 
the instability sets in, and n/2; c) n = 1, at 9 — 0, 0.2587r, at 
the point where the instability occurs, and n/2, for a corner 
junction of ci^2_y2-(-is-wave and s-wave superconductors, with 
length L — lOAj, and zero overlap external current 7°" = 0. 
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FIG. 5. The lowest eigenvalue Ai of the linearized eigen- 
value problem as a function of angle Q, for the n = 0, —1, 1 
solutions. In the range where d is close to zero, the eigenval- 
ues for both n = 0, and —1 are positive and correspond to 
stable solutions. 
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FIG. 6. Overlap critical current 7°" per unit length ver- 
sus the angle Q for a corner junction of d^i_yi ^r is-wave and 
s-wave superconductors, with length L = lOAj, for the vor- 
tex solutions n = 0, —1, 1 that exist spontaneously in the 
junction. 



tices. Let us consider the situation where the junction 
contains a solution in the mode n = 0, at = 0, when 
the net current is maximum. The spatial variation of 
(j) is described by a fractional vortex which is displaced 
around the value (fi = tt, from the corresponding distri- 
bution at zero current which is around 7r/2 (see Fig. |^). 
The current density distribution as seen in Fig. |^a (sohd 
line) at the maximum current is flat above unit with a 
small variation around the junction center giving rise to 
the large value on the net current, seen in Fig. ||. Also at 
= 7r/4 the flat phase distribution corresponding to the 
n = solution at zero current is displaced towards the 
value (p = TT when applying an external current. The cor- 
responding current distribution seen in Fig. (dotted 
line) is straight line and the net current is small for this 
orientation. For the n = —I solution at the point where 
the instability sets in i.e. 9 = 0.2427r, the current density 
distribution is symmetric around zero as seen in Fig. ^ 
(dotted line) and carries zero net current at this point. 
Thus the instability occurs just before the angle where 
a full antifluxon enters the junction. A slightly different 
situation occurs in the magnetic interference pattern of a 
pure s-wave superconductor junction [ p^ where, the net 
current is zero at the magnetic field where a full fluxon 
or antiluxon enters the junction, in the no flux 0-mode. 
At the point = tt/2, of the n = — 1-mode the junction 
contains more than one fluxon causing the characteristic 
oscillations in the current density around the junction 
center as seen in Fig. ^ (dashed line). This reduces the 
critical current for this orientation. 

For the dx2-y2 + id^y pairing symmetry state, we plot 
in Fig. ^) the flux content for the n = 0, —1, 1, versus 
the angle 9. Note the half integer or multiplies value of 
$ at close to or 7r/2. For this grain orientation the 
magnetic flux is only sensitive to the real part of the order 
parameter, which has a sign change but does not break 
time-reversal symmetry. In the d,j.2_y2 -|-is-wave state the 
order parameter is complex for all junction orientations 
and breaks the time-reversal symmetry. Close to or 
7r/2 the flux is fractional. The flux quantization at = 
can be used to discriminate between these states. 

In Fig. ||b) we plot the critical current per unit length 
evolution with the grain angle 9 in the d^2_y2 + idxy-'wave 
state. Close to 6* = we see that the I°" for the n = 0, — 1 
solutions, coincide. This happens also a,t 9 — tt/2 for the 
n = 0, 1 solutions. In these orientations the order pa- 
rameter becomes pure real and does not break the time- 
reversal symmetry. As a result the critical current at 
these angles is the same as in a junction with pure d- 
wave symmetry. At 9 = n/A the order parameter is pure 
imaginary and has the same magnitude for both pairing 
states. As a consequence for 9 — 7r/4, the critical cur- 
rents for both junctions coincide. Also the unstable part 
of the n — 1 branch, in the Ic vs 9 is almost the same 
for the two symmetry states, due to the small difference 
in the flux, compared with the large flux content of the 
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FIG. 7. The current density distribution J{x) of the vortex 
solutions a) n = 0, at ^ = 0, 7r/4, n/2, h) n = —1, ai 9 = 0, 
0.2427r, where the instability sets in, and 7r/2, c) n = 1, at 
— , 0.2587r, at the point where the instability occurs, and 
7r/2, for a corner junction of d^i_y2 + is- wave and s-wave su- 
perconductors, with length L = lOAj, and maximum external 
overlap current /°". 
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FIG. 8. a) The spontaneous magnetic flux $ as a function 
of the angle Q, for the various vortex states, n = 0, — 1, 1, for a 
corner junction between a d^i_yi + idxy-wave and an s-wave 
superconductor, with length L — lOAj. The flux for ^ = 
is integer multiply of $0/2. b) The corresponding critical 
current /°" per unit length. 



solutions in this region. 



B. Magnitude of the secondary order parameter 

In the above calculations the magnitude of the sec- 
ondary order parameter is small compared to the dom- 
inant (i.e. 7120 = O.lriio). However the maximum frac- 
tion of the secondary component that has been observed 
in phase coherent experiments employing different ma- 
terials, geometries, and techniques is up to 25% of the 
dominant This triggered our interest to study the 
magnetic flux and also the critical currents as a function 
of the strength (us) of the secondary order parameter, 
where the magnitude of the dominant order parameter 
Ud is also varied in a way that Us+Ud = 1. When = 
only the d^2_y2-wave order parameter is present, while 
when ris = 1 only the s-wave order parameter appears. 
This situation can be realized for example near the sur- 
face where the dx2_y2-waYe order parameter is suppressed 
and the s-wave order parameter is enhanced. The result 
is presented in Fig. ^) and ||b) for the d^r^-y^ -I- is- 
wave case at 6 = 0. We see that when the secondary 
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FIG. 9. a) The spontaneous magnetic flux <E> and b) 
the critical current 7°" per unit length versus the strength 
Ua of the secondary s-wave component for a corner junc- 
tion of d^2_y2 + is- wave and s-wave superconductors, with 
length L — lOAj, for the vortex solutions n — 0, —1, 1 that 
exist spontaneously in the junction. The magnitude of 
the d^2_y2-wa.ve order parameter is given by the relation 
Us +nd = 1. 

component is absent (i.e. = 0) the picture of the 
drc2 _y2-wave state is reproduced. The same picture also 
holds for the d^2_y2 + idxy-'^s^'ve state at 9 = 0, since the 
order parameter for the dx2-y2 +idxy-waYe state at = 
is real not breaking the time-reversal symmetry. So for 
6^0 the magnetic flux and the critical current for the 
dx2_y2 +idxy-"^ave state would not show any change with 
the variation of the secondary order parameter d^y- As 
Us is increasing the modes n = and n = —1 are no more 
degenerate, in the sense that their flux deviates from the 
value $o/2 and — $o/2 respectively and also their critical 
currents are no longer equal. The mode n = has larger 
critical current because it has smaller flux content in ab- 
solute value. For values of close to unity, the different 
modes contain integer magnetic flux, as in the junction 
between s-wave superconductors, and also their critical 
currents have the same values as in the perfect junction 
problem. The conclusion is that the larger the secondary 
component is in a sample the easier is to be detected in 
a flux measurement experiment. 



C. Magnetic field 

We now examine the influence of the magnetic field on 
the spontaneous vortices for broken time reversal sym- 
metry pairing states. In Fig. 10 we plot the magnetic 
flux at zero current versus the magnetic field H for the 
dx2^y2 + idxy-'wave state at 6 — 0. In the pure s-wave 
superconductor junction there is no overlap between dif- 
ferent modes in the magnetic flux, and each mode has 
magnetic flux which is more than n^Q and less than 
{n + 1)^0- In this problem due to spontaneous magne- 
tization the range of the modes is different and in some 
cases overlapping, and the labeling is with a single in- 
dex n, corresponding to the pure s-wave superconductor 
junction (n,n-|-l) mode p^ . Moreover the range in mag- 
netic flux of each mode is displaced compared to the pure 
s-wave superconductor junction problem by an amount 
which corresponds to the intrinsic flux. Also we have 
the existence of stable vortex states i.e. n = 0,-1, to- 
gether with the unstable ones i.e. n = 1, —2 in a large 
interval of the magnetic field, which is almost the same. 
The n = —2 mode extends to zero magnetic field, and 
the reason we didn't examined this mode in Sec. IV is 
because the stability analysis shows negative eigenvalues 
for all the range of junction orientations, at H = 0. In 
the long s-wave junction the extremum of the mode (0, 1) 
in H is the critical field for one fiuxon (antifluxon) pen- 
etration from the edges, [denoted by Her {Hd), for the 
right (left) edge], and is equal to 2(— 2). The solution for 
the phase at these extremum values of the field becomes 
unstable because the value of the phase at the junction 
edges reaches a critical value. In the problem of a junc- 
tion with some spontaneous fiux, we consider here, the 
range of the corresponding mode in H is significantly 
broadened and also the instability at the boundaries sets 
in due to different reasons. In particular the instability 
occurs due to the interaction of the flux entering from 
the junction edges, when the magnetic field reaches the 
critical value HcriHd), with the spontaneous flux at the 
center. Similar features are encountered in the problem 
of flux pinning from a macroscopic defect in a conven- 
tional s-wave junction. ||l7|] 

We now examine the magnetic-interference pattern for 
the two symmetries where the bias current enters in the 
overlap geometry. In the dx2_y2 -\- idxy-yfave case, where 
9 = 0, this pattern has a symmetric form as we can see 
from Fig. |ll](a) . This is because this result is only sensi- 
tive to the real part of the order parameter, which has a 
sign change but does not break time-reversal symmetry. 
For the angle 9 = 0.5 where the order parameter has a fl- 
nite imaginary part and breaks the time-reversal symme- 
try this pattern becomes asymmetric and the "dip" ap- 
pears to a value of flux slightly different than zero. Note 
that the asymmetry refers mainly to the modes n = 0, 
and n = — 1. The other modes are not influenced much 
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FIG. 10. Magnetic flux $/<l>o at zero external cur- 
rent versus the magnetic field H for a corner junction of 
d^i_yi + id^^y-wsLve and s-wave superconductors, with length 
L = lOAj, for angle 9 — 0° . Hci{Hcr) denotes the critical val- 
ues of the magnetic field where the mode n = 0, terminates. 



due to their higher flux content. Also the critical current 
is suppressed compared to the case where = as can 
be seen in Fig. |l^(b), due to a drop in Jc. 

In the dx2-y2 + is-wave symmetry, in the limit where 
9 —> 0, the order parameter is complex and the pattern 
is asymmetric as can be seen in Fig. ^l]c, for the angle 
9 = 0. This is in agreement with our previous work for 
the inline current input for a junction with (ij;2_,y2 + is 
symmetry jl^ . There it was found that the pattern is 
asymmetric for lengths as long as i = lOAj. For angles 
close to 7r/4, the magnetic interference pattern is similar 
with the dx2-y2 -\- idxy-state. This is because the sin(20) 
dependence of the d^y component is almost unity. This 



is seen in Fig. Hi where we present the variation of the 
critical current per unit length versus the enclosed flux 
for = 0.5, and the symmetry state is d^^-yz + is. 

In the short junction limit L < Xj the same argument 
can be applied without any explicit reference to fractional 
vortex and antivortex solutions. However as we found in 
our previous work [^ , both n = and n = — 1 (there 
/«a, fa) exist, with reduced flux content, in this limit as 
a continuation of the corresponding solutions in the large 
junction limit. In this case the external applied magnetic 
field becomes equal to the self field, and the maximum 
current can be calculated analytically |p^. 



sin(7r$/2<I>o) cos[7r$/2$o + {<Pc2 - (/>ci)/2] 



7r$/2$o 



(11) 



As we see a.t 6 — for the dx2_y2 +idxy-'weive case, the re- 
lation 4>c2—4'ci = ni: holds, and the magnetic interference 
pattern becomes symmetric, while for the d^2_y2-\-is, this 
difference is a fraction of vr and the pattern is asymmetric. 
However as we increase the junction length, we expect 
this symmetric pattern for the d-wave order parameter 
to be continued. This symmetry in the large junction 
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FIG. 11. (a) Overlap critical current per unit length 
versus the magnetic flux "1> in units of <l>o, for a corner junc- 
tion of d^2_y2 +idxy-'wave and s-wave superconductors, with 
length L = lOAj, for angle ^ = 0°. (b) The same as in a) but 
for 6 = 0.5. (c) The same as in a) but for d^2_y2 + is- wave 
and s-wave superconductors for angle 6 = 0°. (d) The same 
as in (c) but for 6 = 0.5. 



limit, is described more effectively by the assumption of 
the n = 0, — 1 solutions which give a symmetric magnetic 
interference pattern as we presented. Also the n — —1 
solution extends to values for the magnetic flux, where 
the n = solution is absent. Eliminating one of them 
will break the symmetry of the diagram. 

VI. EXPERIMENTAL RELEVANCE 

The symmetric pattern with a minimum at zero ap- 
plied field observed in corner junction experiments be- 
tween YBCO and Pb at 9 = has been interpreted as an 
indication of dx2_y2-wa.ve symmetry. |l9| , ^ This result 
refers to short junctions where the junction size is much 
smaller than the Josephson penetration depth. However 
as we found here these experimental data are also consis- 
tent with an order parameter with dg.2_y2 -\- idxy pairing 
symmetry at 9 = 0. 

Also the critical current Ic versus the magnetic flux <i> 
of a SQUID, consisting of two planar Josephson junctions 
on the faces of YBCO superconducting crystal, connected 
by a loop of a second superconductor, for 9 — ot 9 = 
7r/2 is found shifted by $ = 0.5<i>o and has a minimum at 
$ = (instead of a maximum as in a SQUID involving 
conventional s-wave superconductors or the edge SQUID 
in which both junctions are on the same crystal face) but 
is still symmetric. This result has been attributed to an 
order parameter with dx2_y2-waye symmetry. However 
the theoretical analysis done by Beasley et al. shows 
that it is also consistent with an order parameter with 
dx2_y2 -\- idxy-'paiv'mg symmetry at = 0. 

In both cases of SQUID and corner junction the sym- 
metric pattern observed at 9 = Q rules out the dx2-y2-\-is- 
wave pairing state where the order parameter is complex 
everywhere resulting in an asymmetric Ic versus $ pat- 
tern for all angles 9. However the small asymmetry (less 
than 2%) observed at 6* = in some experiments can 
be attributed to various complicating factors e.g. fluxon 
trapping as will be discussed latter in this section. 

The experiment proposed here to resolve ambiguity be- 
tween dx2^y2 + idxy and dx2^y2 at = 0, is to execute 
the same experiments using SQUID or corner junction 
at an angle between sample faces 9 between and 7r/2. 
Our theory predicts symmetric (asymmetric) pattern for 
the dx2-y2-'wave {dx2-y2 + idxy)-'^ave pairing state for 
the corner junction case. This kind of experiment has al- 
ready been done in the case of SQUID geometry. |^2| The 
tunneling directions are defined lithographically and pat- 
terned by ion milling of a c-axis oriented film. A YBCO 
thin film is patterned into a circle with a series of Nb-Au- 
YBCO edge junctions at orientations spaces every 7.5°. 
The measurement of the Ic vs 0, which probes mainly 
the magnitude of the order parameter has an angular 
anisotropy, indicating an anisotropic order parameter. 
Also the execution of this experiments is not easy due 



to the difficulty in cleaning, polishing a crystal at angle 
9, between and 7r/2. 

Also in an experiment analogous to the corner junction 
Miller, Ying et. al. used frustrated thin-film tricrys- 
tal samples to probe the pairing symmetry of YBCO. 
They found a minimum in the Ic vs the externally ap- 
plied fiux $e diagram at <I>e = in the short junction 
limit and a maximum at $e = for a wide junction 
where the junction length is much larger than the Aj. 
However for a wide junction the correct quantity to be 
compared should be the total flux $ which involves con- 
tribution both from the externally applied flux and the 
intrinsic flux. Also in the tricrystal magnetometry ex- 
periments on half-fiux quantum Josephson vorticies one 
can only observe spontaneous magnetization of $o/2 in 
a frustrated geometry only in the large junction length 
limit H. 

There is a number of complicating factors in the inter- 
pretation of the experiments involving corner junctions 
that could lead to an asymmetric (Ic vs $) pattern even 
for 9 = 0. These are the asymmetry of the junction 
(meaning that the critical current of the two junction 
faces are not equal). This will only cause the dip to be 
shallower and will maintain the symmetry of the Ic vs $ 
diagram. Also these experiments are infiuenced by the 
sample geometry and the effect of fiux trapping i.e. there 
can be vortices trapped between the planes of the cuprate 
superconductors that could affect the Ic vs $ diagram. 
In the corner junction case, it creates an asymmetry in 
the fiux modulation curves. However these fiux trapping 
effects are not sufficiently large to change the qualitative 
interpretation of these experiments. 

VII. CONCLUSIONS 

We studied numerically the possible spontaneous vor- 
tex states that may exist in a corner junction between 
a superconductor with time reversal symmetry broken, 
(i.e. dx2^y2 -t- idxy or dx2-y2 + is), and an s-wave super- 
conductor, in the long junction limit. We studied sepa- 
rately three parameters which can be used to modulate 
the spontaneous fiux. These are the magnetic field H, 
the interface orientation 9, and the magnitude of the sub- 
dominant order parameter Us. We pointed out the dif- 
ferences between time reversal broken states under these 
modulation parameters. 

We found that in flux modulation experiments in- 
volving superconductors with some spontaneous fiux the 
range in magnetic fiux of each mode is displaced com- 
pared to the case of a pure s-wave superconductor junc- 
tion by an amount which corresponds to the intrinsic fiux. 
In particular when the magnetic field H is considered as 
the modulation parameter, the range in H of the lower 
fiuxon modes is significantly broadened compared to the 
s-wave case, and the instability at the boundary values of 
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the field sets in due to the interaction of the flux entering 
from the junction edges with the intrinsic flux. In any 
case, for each value of the parameter which changes the 
flux, the modes are separated by a single flux quantum. 

We also derived some simple arguments to discrimi- 
nate between the different pairing states that break the 
time reversal symmetry. For the d^2_y2 + id^^y-wave pair- 
ing state, the junction orientation where 9 = i.e. the 
lobes of the dominant dx^^y2-wave order parameter are 
at right angles for the corner junction, give flux quanti- 
zation condition $ = n<i>o/2 as in the (i3,2_j,2-wave state, 
which is different from the corresponding flux quantiza- 
tion for the dx2-y2 -t- is-wave pairing state, at = 0, 
which is <I> = (n/2 -|- /)<i>o, where / is a small quantity. 
These different conditions provide a way experimentally 
to distinguish between time reversal broken symmetry 
states. Note that since the magnitude of the secondary 
order parameter is small compared to the dominant, the 
detection of time reversal broken states requires a very 
precise measurement of the spontaneous magnetic flux. 

Also we showed that the magnetic interference pattern 
at = is symmetric (asymmetric) for the d^2_y-2 -\-idxy 
{dx2-y2 + is), and this also can be used to probe which 
symmetry the order parameter has, at least where the 
junctions are formed. We expect our findings, for the 
magnetic field dependence of the critical current, to hold 
even in the short junction limit, where the most experi- 
ments on corner junctions have been performed PJl^. 
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